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Abstract 



We study a reflectionless "PT-symmetric quantum system described by the pair of 
complexified Scarf II potentials mutually displaced in the half of their pure imaginary 
period. Analyzing the rich set of intertwining discrete symmetries of the pair, we find 
an exotic supersymmetric structure based on three matrix differential operators that 
CN " encode all the properties of the system, including its reflectionless (finite-gap) nature. 

The structure we revealed particularly sheds new light on the splitting of the discrete 



If} . states into two families, related to the bound and resonance states in Hermitian Scarf 

II counterpart systems, on which two different series of irreducible representations of 
■ sZ(2, C) are realized. 

^ \ Keywords: Supersymmetric quantum mechanics; VT -symmetry; non-Hermitian Hamil- 

tonians; self-isospectrality; finite-gap systems; non-linear super symmetry. 

^ : 1 Introduction 

. 

Fourteen years ago, Bender and Boettcher discovered a huge and remarkable class of non- 
Hermitian quantum Hamiltonians which exhibit an entirely real spectrum |T] . One of the key 
points of the observation was the requirement of a VT-symmetry generated by the product 
of the parity, V, and time, T, inversion operators, 

VxV = —x, TxT = x, TiT=—i, (1.1) 

that substitutes the usual quantum mechanical property of Hermiticity. This condition, 
however, is necessary but not sufficient for the reality of the spectrum; it is also required 
that eigenfunctions of a PT-symmetric Hamiltonian must be simultaneously eigenfunctions 
of the VT operator. In this case we say that the VT symmetry is unbroken, otherwise the 
eigenvalues are not real, a part or all of them appear in complex conjugate pairs, and the 
VT symmetry is broken [H [21 [3] . 
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Nowadays, quantum mechanics with non-Hermitian Hamiltonians transformed into an 
independent line of research where, specifically, the notion of the PT-symmetry was gener- 
alized into the condition of pseudo-Hermiticity [I]. Non-Hermitian Hamiltonians appear in 
physics in diverse areas including quantum optics, cosmology, atomic and condensed matter 
physics, magnetohydrodynamics, among others. For a good review of the developments in 
the area and applications, see [H |6] and references therein. 

On the other hand, there is a wide class of modern techniques and methods which prove 
their effectiveness in the study of quantum mechanical systems. One of them is supersym- 
metric quantum mechanics (SUSYQM), introduced initially by Witten as a toy model to 
study the spontaneous supersymmetry breaking [7j. Over the past few decades, SUSYQM 
transformed into a powerful tool in quantum physics [8] which turns out to be useful, for 
example, in the spectral analysis as well as in searching for new solvable systems. SUSYQM 
in its usual form is based on a Darboux transformation that relates two Hamiltonians by 
means of an intertwining, linear differential operator, and leads, as a result, to a complete 
or almost complete isospectrality of both systems [9] . 

The concept of Darboux transformations and SUSYQM was generalized in different as- 
pects that lead to the discovery of the classes of quantum systems that reveal non-linear 
[Ha EB E2], bosonized PU EH El US] and self-isospectral \JM El HB3 Q3U ED] supersym- 
metries. Particularly, a certain class of potentials was found in which all the mentioned 
specific types of supersymmetries were brought to light in a form of a peculiar structure 
that was coined in [20] as u tri-supersymmetry n . Such a structure was shown to underlie 
special properties of some physical quantum systems [211 E2J E3j EH |25] . Among the prin- 
cipal characteristics of tri-supersymmetry, which will be described below, is the existence 
of three integrals of motion in the form of supercharges which encode the main properties 
of the corresponding systems. An example of the systems that reveal a tri-supersymmetric 
structure is provided by the Hermitian finite-gap periodic potentials [20], which in the limit 
when their real period tends to infinity are known as reflectionless potentials. 

The idea of supersymmetry in some of its versions was adapted in the context of non- 
Hermitian Hamiltonians [23 ESI [221 E01 EU E21 E3] - Although the non-linear and self- 
isospectral supersymmetries were studied before in the systems with non-Hermitian Hamil- 
tonians [31] , the presence in them of a supersymmetric structure that would unify the men- 
tioned types of supersymmetry remains to be unknown. One can wonder therefore if there 
exist non-Hermitian potentials of finite-gap nature that display the properties to be simi- 
lar to those of the Hermitian counterpart, and what are the peculiarities of the associated 
supersymmetric structure. 

The purpose of the present article is to report the observation of an exotic tri- 
supersymmetric structure in a broad class of PT-symmetric extended finite-gap systems. 
This allows us, on the one hand, to understand and explain the main features of their 
spectra by analyzing the properties of the supercharges; on the other hand, we explicitly 
trace out the differences that appear in comparison with the Hermitian case, at the level of 
potentials as well as in the generic properties of tri-supersymmetry. The importance of a 
hidden, pure imaginary period in non-periodic on a real line finite-gap systems is clarified, 
particularly, in the light of tri-supersymmetric structure. We achieve all this by considering 
the pairs of PT-symmetric complexified Scarf II potentials mutually displaced in the half of 
their unique imaginary period. For special values of the parameters these potentials become 
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perfectly transparent, i.e. reflectionless. It is exactly for such special parameter values the 
tri-supersymmetric structure arises and allows us to clarify its nontrivial interplay with var- 
ious discrete symmetries of the potentials. The supersymmetric structure we reveal gives us 
a new insight on the splitting of discrete states in such a class of PT-symmetric quantum 
systems into two distinct families, on which two different representations of the sl(2, C) alge- 
bra are realized, the fact that was established initially by Bagchi and Quesne by using group 
theoretical methods [35] . Note here that the complexified Scarf II potential was extensively 
studied in the literature in the context of PT-symmetric and pseudo-Hermitian Hamilto- 
nians; the updated summary on these investigations can be found in ref. [36] • Recently, 
this potential also attracted the attention in different areas of physics such as quantum field 
theory in curved spacetimes [37], soliton theory in nonlinear integrable systems [38] and also 
the physics of optical solitons [39] . 

The plan of the article is as follows. In sections [2] and [3] the basic properties of the pair of 
complexified mutually conjugated Scarf II potentials are reviewed. Specifically, in the next 
section we study various discrete symmetries of the potentials, describe spectral properties 
of the systems in dependence on the parameter values, and describe the relation with the 
other known potentials. In section [3] we first present the explicit expressions for the two 
families of the singlet states in the spectrum, and then discuss the continuous spectrum and 
its relation with that of the free particle by means of Darboux-Crum transformations that 
underlie the non-linear supersymmetry. The tri-supersymmetric structure is described in 
section [4j Section [5] provides two concrete nontrivial examples of the systems with two and 
three bound states to illustrate the general results. In section |6] we present the discussion 
and concluding remarks. 



2 Discrete symmetries and relations 

Consider a pair of complexified Scarf II potentials 

+ , s I 2 + m(m + 1) , sinhx 

V?J X ) = b -±« 2m + l — -s- . (2.1) 

cosh x cosh x 

We suppose that / and m are real parameters, and x G R. Potentials (12. ip are then free of 
singularities on the real line, and their real and imaginary parts vanish for x — > ±oo, see 
Fig. [H Pair (12 .11) can be transformed into a pair of original, real Scarf II potentials [40] by 
a substitution I — >• il, see Fig. [2j Both the original and complexified potentials have a pure 
imaginary period 2iir. 

Parity inversion, V, and time reversal, T, operators intertwine the completely isospectral 
potentials (12. ip . 

VV,% = V^V, TV t % = V^T. (2.2) 

As V^i m = y^_ m __ x = V^, another pair of intertwiners is provided by operators IZi and lZ m , 
which act on the parameters / and m, TZi : (/, m) — > (— /, m), lZ m : (/, m) — )■ (/, — m — 1), 

KK m = V^Tlu H m V± m = V^JZ m . (2.3) 
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Figure 1: Plot of real and imaginary parts of the complexified Scarf II potential Vf^, which 
here as well as in a generic case of (12. ip are even and odd functions, respectively. 
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Figure 2: Examples of original (real) Scarf II potentials. 



Yet another pair of intertwining operators is given by pure imaginary translations for the 
half of the period 2?7r, T + : x — > x + iir and T_ : x — > x — in, 

T + V t % = Vf m T + , T_V t % = Vf m T_ . (2.4) 

The product of any two of the listed intertwining operators, except of T + = exp^^) and 
T_ = exp(— 7T^), T + T_ = 1, is a nontrivial discrete symmetry of each of potentials (12. ip . 
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Particularly, we find that each potential of the pair is 'PT-symmetric, 



[PT, Vfj = . (2.5) 

The potentials satisfy also the relation 

V t % = V^-m-l , (2-6) 

which is produced by a composition lZ(IZ m of the intertwining generators (12. 31) . 
We also have a symmetry relation 



^ = (2 - 7) 



that will play a key role in the analysis below. The composition of (12. 6 1) and ( 12. 7p produces 
yet another symmetry 

V Cm = V -m- M ■ ( 2 - 8 ) 



2 ' 2 



The nature of relations H 2 . 3 f) and ( I2.6l) -( l2~8l) takes a somewhat more transparent form if to 
redefine the parameters [41] : (/, m) — > (a + , a_), a-t = m±/ + |. Then the coefficients in (12.11) 
are transformed into / 2 +m(m+l) = |[a++a 2 — |] and/(2m+l) = \{a 2 + — o?_ ), and symmetry 
(12. 7p corresponds to a reflection V a _ : (a + ,a_) — > (a+, — a_) in the plane of a-parameters. 
Intertwining relations ( 12. 3 p are given by the products of reflection V a _ and a-rotations for 
±tt/2, TZ +7t/2 : (a+, a ) -> (a_, -a+) and 11^/2 ■ («+, «-) -> (-«_, a+), ^ = V a _K +7T/2 , 
Tim = V a _7Z- 7T /2, while symmetries (12. 6p and (12.81) correspond, respectively, to a 7r-rotation 
(a + , a_) — > (— a+, — a_) and to a reflection in a+, V a+ : (a + , a_) — > (— a + , a_). 

It will be more convenient for us, however, to work in terms of the parameters / and m. 
Relations ( I2.6p -( 12T81) allow us to interchange integer values with half-integer ones as well as 
positive with negative values. By this reason, without loss of generality we can suppose that 
/ and m are non-negative integers; sometimes, however, negative and half-integer values will 
be important as well. 

The described discrete relations and symmetries form the base for a rich exotic super- 
symmetric structure that we will reveal in the extended system = diag(Hf m , Hf m ), 

H tm = -^ + V l %, (2.9) 

with pseudo-Hermitian Hamiltonians 

Ht=V{H±JV, (2.10) 

where we have taken into account the relation P" 1 = V . 

The spectra of H lm display different features in dependence on the values of the param- 
eters / and m. The eigenf unctions of Hamiltonians (12. 9p may or may not be simultaneous 
eigenstates of the VT operator. The situation is known as unbroken VT symmetry in the 



1 A pair of complex potentials ()2.1|) with m = and corresponding extended Hamiltonian T-Li^ emerged 
recently in |37j under investigation of quantum field theory in dc Sitter space. 
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former case, or spontaneously broken in the latter case. As was shown in [TJ [2j [3], the spec- 
trum of a VT symmetric Hamiltonian is real when its eigenf unctions are simultaneously the 
eigenstates of the VT operator. In the case of the H^L , this situation holds when 

I 2 + (m + 1/2) 2 > \2l(m+ 1/2) | , (2.11) 

that is always valid for real / and m [32]. For real / and m, the spectra of H^ m and Hf 
have a finite number of bound states of non-degenerate energies, while the probability flux is 
conserved in a scattering sector. It was pointed out in [32] that an imaginary displacement 
in the spatial coordinate x — > x + i5, where iJ G K, breaks in general the conservation 
of probability, i.e. |T| 2 + |i?| 2 ^ 1. There is, however, a case that shows several special 
properties for this class of shifted potentials. When I and m take simultaneously non- 
negative integer values, the reflection coefficient vanishes, \R\ = 0, and the probability is 
conserved independently of any displacement 5 ^ ±tt/2, see below. As a result potentials 
(12. ip belong to the class of reflectionless potentials [33]- In this case both potentials have the 
same spectrum with m + l + 1 singlet states, where m + / of them are bound states of negative 
energies, and there is a singlet state of zero energy at the bottom of the continuous part of 
the spectrum. We discuss the spectral characteristics of (12 .9p in more details in Section [3j 

Within a framework of PT-symmetric quantum mechanics [5], or in a more general 
framework of the systems with pseudo-Hermitian Hamiltonians [6], it was shown that the 
VT integral can be supplied with yet another nontrivial (nonlocal) integral of motion that 
has a nature of the charge conjugation operator, C [3S]. This allows finally to define a 
positive definite scalar product to extend a probabilistic interpretation for the case of VT- 
symmetric systems. Here, we just have in mind this general picture when discuss the bound 
and scattering states by referring to original papers on the subject. 

The double degeneracy in the spectrum for the scattering sector together with the reflec- 
tionless property and the finite number of singlet states indicate that each system Hf m and 
Hf m possesses a hidden, bosonized non-linear supersymmetry [ITJ [13] as this happens for 
Hermitian reflectionless Hamiltonians [T^| [20] . In fact for integer values of I and m, poten- 
tials (12. ip are solutions of the KdV hierarchy: they satisfy the stationary s-KdV n , n = l + m, 
non-linear equations [36] . This means that there should exist non-trivial integrals of motion 
^2n+i an d ^2n+i f° r Hamiltonians (12 .9p in the form of differential operators of order 2n + 1. 
These integrals should satisfy relations 

[A£*.i, H+J = [A 2 - n+1) H-J = 0, (A± +1 ) 2 = P(fljfJ, n = l + m, (2.12) 

where P(H^ m ) is a polynomial of order In + 1, i.e. A2 n+1 and H^ m , as well as A^" n+1 and 
H[~ m , should compose a Lax pair [33 H8] . 

There are particular cases for which H^ m are completely isospectral to the reflectionless 
Poschl- Teller Hamiltonians I49I 2 I 



x 



2 The free particle corresponds to A = and A = — 1, and can be considered as a zero-gap case of the 
family of finite-gap, reflectionless Poschl- Teller systems ()2.13|) [14l [50]. 
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Namely, the Poschl- Teller and the complexified Scarf II potentials are related by a transfor- 
mation of an imaginary displacement and a rescaling of the variable, 



m = I — 1, A = m + 1, 
iff T (x ± ztt/4) = 4H± m (2x) , for <{ and (2.14) 

m = I, X = m . 

Another interesting relation is with a generalized, singular Poschl- Teller potential [8] 

= t 3 +m ( ?+ i ) + ;(2m + l) ^ 

sinn x sinn x 

The singularity can be removed by complex shifting [35] • For special values of such a shifting 
we get the PT-symmetric pair of potentials (12. ip . 

V l % PT (x±z7r/2) = V l %(x). (2.16) 

It is worth to note that the potentials ( 12.11) and corresponding Hamiltonians admit a 
representation that generalizes Eq. (I2.14p . 



H+ - * r(r + 1) S(S + 1) \ (217) 

l > m 4 1 de cosh 2 e cosh 2 (e+^^ 1 ' 




1 _f_ s(s + l) r(r + l) 
l ' m 4 ^ de cosh 2 e cosh 2 (e 1 

where £ = | + i|, and 

r = / + m, s = I — m — 1 . (2-19) 

Representations (I2.17P and (I2.18P correspond, on the one hand, to the relation (12 .4p between 
potentials (12. ip generated by the shift in the half of their imaginary period, 

H^) = H tm{i + ^) ■ (2-20) 

On the other hand, the generators of intertwining relations (12. 3p correspond here to TZi : 
(r, s) — > (—s — 1, — r — 1) and lZ m : (r, s) — > (s, r). 

We will return to ( I2.20p in the discussion of the supersymmetric structure, but here 
we note that in the context of periodic (elliptic) finite-gap potentials, a so called tri- 
supersymmetry appears when the superpartner potentials are the associated Lame potentials 
shifted mutually in the half of the real period [20] • As we will show below, the extended 
system T-L^ m constructed from (12.91) also possesses a tri-supersymmetry in which the imagi- 
nary period, the VT- symmetry, the discrete symmetries ( 12. 6 p and (12. 8p . and the non-linear 
supersymmetry together play a fundamental role to form altogether a unified structure. 

As a final comment on relation of ( 12. ip with periodic finite-gap potentials, we note 
that there is a generalization of the family of the associated Lame potentials known as the 
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Darboux-Treibich-Verdier potentials [51]. In terms of the (double periodic) Jacobi elliptic 
functions they read 



V DTV = n 1 (n 1 + l)k 2 sn 2 x + n 2 (n 2 + l)k 2 sn 2 (x + iK') 

+ n 3 {n 3 + l)k 2 sn 2 {x + K + iK')+n 4 {n 4 + l)k 2 sn 2 {x + K) (2.21) 

= nAni + l)fc 2 sn 2 a; + n 2 (n 2 + 1) — 5 — h n 3 (n 3 + 1) — - — h n 4 (n 4 + 1)— — . 

sn 2 x cn 2 x dn x 

Here < k < 1 is the modular parameter, and ( I2.21j) has a real, 2K, and an imaginary, 
2iK', periods; K = K(k) is the elliptic complete integral of the first kind and K' = K(k'), 
k' = \/l — k 2 [53]. The finite-gap nature of (I2.2ip appears when parameters rii take integer 
values. Particularly, when n 2 = n 3 = 0, potential (I2.2ip reduces to the finite-gap associated 
Lame potential [20j |52]. When the modular parameter takes the limit k — > 1, the real period 
tends to infinity, 2K — > 00, while 2iK' — > in, and the potential transforms into 

V™ — + + - n2( " 2 + 1 \ + const , (2.22) 

cosh a; cosh (x + z^) 

that has the form of potentials (12.171) and (12.181) . 

In the next section we will study the states of the PT-symmetric systems Hjj^ in the 
light of the discrete symmetries and relations that we have discussed. 



3 Wavefunctions and differential intertwiners 

The group theoretical methods and supersymmetry are the powerful tools in the study 
of quantum mechanical systems. The Hamiltonians f ]2.9f) provide a good example of the 
systems for which these techniques work effectively, particularly, to analyze the spectrum and 
eigenf unctions. In this direction, using irreducible representations of the s/(2, C) algebra, it 
was found in [35] that the non-degenerate parts of the spectra of potentials (12. ip are described 
by the two sets of eigenfunctions, one of which ig^| 

= sech m x exp fril arctan(sinhx)]P^ m - 1/2 '-^ m - 1/2 (±2 sinhx) , (3.1) 

where P n Q,/3 (x) are the Jacobi polynomials [53]. The corresponding energy levels and the 
values of the parameter n are 

E n ,m = -ijn - n) 2 , n = 0,l,2... <m. (3.2) 

3 It is worth to note that in the early stages of studying the complexified Scarf II potential, just one series 
of the singlet states, (|3.1|) . that comes by analytic continuation of the Hermitian version, was considered in 
the literature [281II2]- Later, the complete set was found by algebraic methods in [35] and reconfirmed in 
refs. [41] [43] [54] . The second, lost set of states, p.4[) . corresponds to resonances in the Hermitian counterpart 
potential. Within the problem for the complexified Scarf II potential this second set can be obtained from 
the counterparts of the singlet states of the Hermitian problem by applying the symmetry transformation of 
the parameters (|3.3j) . 
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The eigenfunctions ^f^m for n < m describe m bound states, while n = m corresponds to 
the singlet zero energy state at the bottom of the continuous spectrum. With the discrete 
symmetry (12.71) of the potentials, Vf^ = V i i: to which corresponds a transformation 
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(l,m) -> (m + i,/-|) , (3.3) 
it is possible to write down the another set of the bound states, 

= sech'- 1/2 x exp [^i(m + 1/2) arctan(sinhx)]P™ +1/2 -''- m - 1/2 -'(±i sinhx) , (3.4) 
with energies 

E n ^ = -{I - n - 1/2) 2 , n = 0,1,2... </- 1/2, (3.5) 

so that (13. ip and (13. 4p represent together the I + m + 1 singlet states of the systems (12. 9p . 

As we will see, the separation of singlet states of each subsystem Hf m and H[~ m into 
the two subsets is reflected by a specific nonlinear supersymmetry of the extended system 
TLi }m . This supersymmetry is related to the (imaginary here) mutual half-period shift of the 
subsystems. To the best of our knowledge a similar kind of supersymmetric structure, that 
we shall discuss in the next section, was discussed till the moment only for finite-gap systems 
with Hermitian Hamiltonians [20]. To explain this structure and its origin, we present below 
some further comments on the properties of the Hamiltonians (I2.9p and their eigenfunctions. 

States (13. ip and (13.41) can also be obtained from a supersymmetry approach, by means of 
the Darboux-Crum transformations (for details we refer to [9]), so that it is possible to link 
the Hamiltonians with different values of / and m between them. Particularly, it is possible 
to relate the free particle, for which H^ = — = H , with the generic case H^ m . Note that 
the free particle system is presented equivalently here also by Hq_ 1 , Hy 2 and H^^ 2 _ 1 , 2 . 
The bound states described above can be computed from the appropriate non-physical states 
of the free particle, whereas the states from the continuous part of the spectrum are obtained 
from the plane wave states of the free system. To illustrate this picture we will show how 
the scattering sector of H* m can be obtained having in mind that all the results for Hf m 
can be reproduced then with the help of symmetries (I2.2p . or by the shift for the half of the 
imaginary period. The construction of the bound states from the non-physical states of the 
free particle will be illustrated in Section 

Let us define the first order differential operators 



A ± = — ± 

l > m dx 



I | tanhx + i ( m H — I sechx 

2 \ 2 



(3.6) 



They generate the intertwining relations 

A,m^,m = ^i-l,mA,m ' A,m^-l,m = ^mA.ra ' (3-7) 

Now, by applying, as in the singlet states case, the discrete symmetry (13.31) to ( 13. 6ft . we 
obtain the operators 

B, m = — ± (mtanhx + z/sechx) . (3.8) 
dx 

Application of symmetry (13. 3p to (13. 7p results then in the intertwining relations 

Bj m H lm = H t m , B l H l x = H t B, m . (3.9) 



9 



The operators Af m are related between themselves by Af m = V{Af m )^T > 1 . Similarly, for 
the operators B^ m we have Bf m = ViBf^V' 1 . Such relations of conjugation underly the 
pseudo-supersymmetry discussed in the literature for non-Hermitian systems [33], particu- 
larly, a PT-symmetric one. 

Coherently with the discrete symmetry (I3.3p . operators (I3.6P and (13. 8 p allow us to fac- 
torize, up to an additive constant term, the same Hamiltonian in two different ways [53] , 

= -Bl m Bt m -m\ (3.10) 

Since A± m as well as Bf m are PT-antisymmetric (PT-odd), {PT, Af m ] = {VT,B± m } = 0, 
the Hamiltonians Hf m are the PT-symmetric (PT-even) operators. 

In the next section we will see further implications of existence of these two related 
types of factorization operators. For the moment, it is worth to note that the action of the 
operators Af m and B^ m on the two-parametric family of Hamiltonians (I2.9P is quite simple: 
while the former act by lowering and raising the parameter /, the latter play the same role 
for m. Using this fact, we are able now to connect, for any I and m, the Hamiltonian H^ m 
with the free particle Hamiltonian H in two different ways by making use of the operators 

2\m = B l m B l m _ x . . . B l2 B ll A lQ A l _ lQ . . . A 2fi A 10 , (3.11) 
2\m = B^ l _ m B^ l _ m+1 ...B_ l _ 1 B_ lfi A_ l+10 A_ l+2fi ...At 1>0 A 0fi . (3.12) 

The differential operator T>f m has here the order m+l, meanwhile P[~ m has the order m+/ + l. 
These operators intertwine Hamiltonian H^ m with the free particle Hamiltonian, 

2>i7m#o = Kn V T, m . Vr m H = H+ m V Tm . (3.13) 

From (I3.13P the inverse intertwining relations are easily obtained by defining T>^ m = 
n^xJ V ~ X and Hm = -V&iJV- 1 , that yields 

Kn H tm = H oV+ m , V+ m H+ m = H V+ m . (3.14) 

Several comments are in order here. First we note that the difference of orders of the 
two intertwining operators (13.111) and of their conjugate ones are related to the fact that the 
system H^ m can be presented alternatively by the equivalent Hamiltonian Hi_ m _ 1 . Another 
point is worth to note is that since operators A^ m and B lm do not commute, and the path 
that connects the points (I, m) and (0, 0) in the parameter plane can be chosen in different 
ways, the corresponding intertwining operators have not a unique form. It is the operators 
r D\ m and T)f m that together with the corresponding pair of Hamiltonians H^ m will form a 
basis for the construction of the tri-supersymmetric structure in the present "PT-symmetric 
case, see Figure[3l 

On the other hand, the fact that two different Darboux-Crum transformations can relate 
two different quantum mechanical systems is known for the case of Hermitian operators and 
was exploited in [20] to reveal a peculiar, tri-supersymmetric structure in some periodic and 
non-periodic finite-gap systems. Particularly, it was shown in [23] that the two non-trivial 



H + — —A" A + 
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Figure 3: The family of complexified reflectionless Scarf II systems may be presented on 
integer or half-integer lattices in the (l,m) parameters plane [in axes a± = m ± I + |, 
the lattices for Hf m have half-integer coordinates]. Four different points correspond to the 
same system, two of which, shown here as an example by the filled circles, are on the 
integer lattice; another two points, shown by the unfilled circles, are on the half-integer 
lattice. The equivalent points on the same lattice are related by the symmetry transformation 
(12.61) . while equivalent points on different lattices are related by the discrete symmetries 
(12. Tl) and ( 12.81) . The filled and unfilled triangles correspond to a system shifted in the 
half of the pure imaginary period; two such mutually displaced systems are related by the 
intertwining discrete transformations (12.21) and (12. 3p . The two filled and two unfilled squares 
correspond to a free particle system. Any two systems represented on the same lattice may be 
related between themselves by differential intertwining operators. Particularly, any nontrivial 
complexified reflectionless Scarf II system may be intertwined with the free particle. Two of 
such (of many possible) " intertwining paths" shown for the system H^i = H^ 2 _ 2 correspond 
to the action of the operators (13.111) and (13. 12[) with I — 2, m = 1. 



Darboux-Crum transformations encode the existence of a Lax pair in reflectionless Poschl- 
Teller systems. As we noted above on Eq. (12.141) . there are cases in which the potentials 
( 12. ip are reduced exactly to the shifted Hermitian Poschl- Teller potentials, so that operators 
( 13. lip and ( 13 . 1 2[) match in those particular cases the corresponding operators in [201 150] . 

One of the direct applications of the constructed intertwining operators is that we can 
use them to map the plane wave states of the free particle, 

H ^ ±k = k 2 */j ±k , t/j ±k = e ±lkx , (3.15) 

into the scattering eigenstates of Hf m , 

HtA = k 2 n„ n k = K m e ±lkx = c^V-^* , (3.16) 
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where k > and C\{k) is some (fc- dependent) constant factor. For k = 0, the action of the 
operator X> z _ m produces the unique singlet state of the continuous spectrum of , and then 
( I3.16P coincides with the eigenfunction ( 13.1 p with n — m, i.e. = ^^ n - O n the other 
hand, X^~ m annihilates the singlet state located at the bottom of the continuos spectrum of 
the free particle, i.e. ci(0) = in (13.16j) . 

This picture of the Darboux-Crum transformations explains, as in the case of Hermitian 
finite-gap systems [20], the reflectionless properties of the Hamiltonians (12. 9p for integer (and 
half-integer) values of / and m. 



4 Tri-supersymmetric structure 

In this section we will show how the PT-symmetry, originated from the intertwining relations 
( 12. 2p . the discrete symmetries behind (12. 3p . the self-isospectrality based on (12 .4p . the Lax 
integrals A 2n+1 , and the non-linear supersymmetry form altogether a peculiar structure. 

To reveal and describe such an unusual extended nonlinear supersymmetric structure, we 
will show first that the extended PT-symmetric Hamiltonian 

U,, m = ( f ^ ) (-) 

for / ^ has three mutually commuting non-trivial basic integrals of motion, anti-diagonal 
Xi t m and y^ m , and diagonal Z^ m = diag (Aj n+1 , A^+i), n = I + m, 

[Xl,m, Hi,™] = 0, [yi t m, Hl,m[ = 0, [Zl,m, Hl,m] = . (4.2) 

These are the matrix non-linear differential operators of the orders \Xi tm \ = 21, \y^ m \ = 2m+l 
and \2>i, m \ = 2(/ + m) + 1, connected between themselves by a factorization relation 

Zi,m = x^ m yi )m = y^mX^m ■ (4.3) 

It is due to these three basic nontrivial integrals the corresponding supersymmetric structure 
is referred to as a tri- supersymmetry. We will see that it reflects coherently the peculiar 
properties of the extended complexified Scarf II system 1-Li jm , including the existence of two 
types of the discrete energy levels in its spectrum. 

The case I — is particular since the extended Hamiltonians 04. ip just reduce to the 
two copies of the Hermitian Poschl- Teller systems with coupling parameter A = m in ( I2.13p . 
i.e. Hq X = Hq X = H^ T . One integral then reduces to the Pauli sigma matrix, X 0tm = o\. 
The remaining two integrals are related in accordance with 04. 3 p as 3^o,m = o~iZ 0)Tn , where 
the diagonal matrix elements of Zo,m coincide, A^" m+1 = A 2m+1 , and generate the hidden 
bosonized nonlinear supersymmetry of the reflectionless Poschl- Teller system H^ T , see j!4j . 

Before passing over to the construction of the nontrivial integrals, we note that the 
extended system 04. ip is formed by the two self-isospectral Hamiltonians displaced mutually 
by the half of their imaginary period. As a result, the degeneracy of the spectrum of Jii 7m 
is twice that of its corresponding diagonal components. As in usual (Hermitian) quantum 
mechanics, by virtue of relations 04. 2p . it is natural to expect that there is a basis where all 
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the eigenstates of the Hamiltonian (14. ip are also the eigenstates of the nontrivial integrals 
of motion A^ m , 34m, and Z^ m . In accordance with this, as we will see, the m + I doublet 
states corresponding to the set of bound states can be presented in the form 



<i= ( ), n = 0,l,...,</-l/2, (4.4) 



< m =( 2h ), n = 0,1,2... <m, (4.5) 



and 

$± = ( 

\ m,l 

with eigenvalues given by (13. 5p and (13.21) . respectively. The scattering states can be written 

as 

U*:J' (4 ' 6) 

with energies E = k 2 , k > 0. The energy levels with E > are then four-fold degenerate, 
while for E = we have $^ = = and, as in the bound states case, the double 
degeneration. Up to a multiplicative constant, $^ coincide with (14 .5p with n = m. 

The intertwining relations (12.31) for the potentials are trivially extended for the Hamilto- 
nians fT2"H . 

^ m = H^JZ h TZ m Ht m = H^ m TZ m , (4.7) 

where the generators IZi and lZ m can be used to construct a discrete symmetry lZ{R m for 
the extended Hamiltonian l-ii. m - For the extended system Hi, m , antidiagonal Pauli matrix 
o\ (as well as 02), produces the same effect of intertwining of the Hamiltonian's components, 
<j\diag (Hf m , Hf rn )ai = diag (H[~ m , Hf m ). Therefore, we can construct the matrix operators 

Til = aiJli, TZm = a{R m , (4.8) 
which are the integrals of motion for our extended system (14.11) . 

[nt,ni, m ] = 0, [n m ,ni, m ] = o. (4.9) 

In the case of TZi, the commutation relation comes from the intertwining relation, which, in 
turn, is based on the equality B._ x m = Hf m . In the previous section we have seen that the 
operators Af m acting on the Hamiltonans H lm , can lower or raise the index I, see Eq. (13. 7p . 
by means of a chain of Darboux transformations. This means that the appropriate product 
of the operators Af m produces exactly the same intertwining effect as the IZi, which changes 
(l,m) for (— I, m). Indeed, this can be achieved by application of the 2/-th order differential 
operators 

^i, m = A_ l+l m A_ l+2tm . . . A^ m . . . A l _ l rn A l m , (4.10) 
-^■i,m = A lim A l _ lrn ...A 0m ...A_ l+2m A_ l+lm , (4-11) 

which satisfy the intertwining relations of the same form as in (14 .7p . 

X Lm H Lm = H -l,m X Lm = H Lm X L m ■ ( 4 -12) 
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The case I = reduces trivially to the operators X^ m = 1. On the other hand, the nontrivial 
analog of the discrete symmetry TZi, TZf — 1, is provided by the matrix differential operator, 

" ( *i X o m ) • (4 ' 13) 

which, by virtue of (I4.12p . is an integral of motion. In addition to (14. 2p . the integral (14.131) 
satisfies a superalgebraic-type relation 

{Xl,m, Xl,m} = 2X* m = 2P X (U l , m ) , (4.14) 

where P% is a polynomial of order I in Hamiltonian 1-t^ m , 

[1-1/2] 

Px{-H Lm ) = \{ (H l>m + {l-8- 1/2) 2 ) 2 . (4.15) 

s=0 

In the context of analogy of X^ m with TZi, Eq. (I4.14p is a generalization of the relation TZf = 1. 
The polynomial P% has the nature of a spectral polynomial, but which only includes the 
energies of one set of bound states (14.41) . From here a remarkable property of Xi m can be 
derived: acting on eigenstates of "H; im , it annihilates one complete set of doublets while the 
doublet states of another set are the eigenvectors with nonzero eigenvalues, 

[1-1/2] 

X Lm ^ = 0, X hm <S>t, m = ±(-1)" II ^™ + ( l ~ s - V2) 2 )< m , (4.16) 

for n = 0, 1, ...,</— 1/2 and n — 0, 1, 2, . . . < m. Therefore, the integral X[ m identifies all 
the states which correspond to resonances of the Hermitian counterparts of the Hamiltonians 
H^ m [13]. The action of the operator X^ m on the scattering states is characterized by the 
property that it does not distinguish the waves coming from the left or from the right, but 
separates the states with distinct values of the upper index, 

[1-1/2] 

= ±(-i)' n (^+(/-s-i/2) 2 )$j fc , 

s=0 
[1-1/2] 

X±?t k = ±(-1)' II (k 2 + (l-s-l/2f)^_ k . (4.17) 

s=0 

We can construct also differential operators of order 2m + 1 , 

Y Cm - B t-m B t-m+l ■ ■ ■ B tfl ■ ■ ■ B tm-l B tmi ( 4 -18) 
Y l,m = B l,m B l,m-l ■ ■ ■ B l,0 ■ ■ ■ B l,-m+l B l,-m ' (4.19) 

which generate the Darboux-Crum transformations similar to the intertwining relations pro- 
duced by lZ m , 

Y l,m B Lm = B L-m-l^Lm = ^LrrXhrn • (4.20) 
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With their help, we find that the matrix differential operator 

*- ^ ( i ) (4 ' 21) 

is the another nontrivial integral of motion for the extended system 1-L^ m . Like the 1Z m 
commutes with the TZi, the nontrivial integrals (I4.2ip and (I4.13P also commute, 

[X hm ,yi, m ] = 0. (4.22) 

The integral 3^z,m generates a relation 

m—1 

{^,m, yi, m } = 2Py(Hi, m ) = 2U hm JJ (ft,, m + (m - r) 2 ) 2 (4.23) 

r=0 

to be of the form similar to (I4.14p . The roots of the spectral polynomial Py(Tii,m) ar e 
complementary to those of the polynomial Px{Hi,m)'- they coincide with the energies (13.21) 
of the eigenstates (14.51) . In correspondence with this property, the second non-trivial integral 
of motion, 3^,m, annihilates the remaining set of discrete eigenstates of H^m, not annihilated 
by the integral Xi yTn , which correspond to the doubly degenerate energy levels, while the zero 
modes of the latter integral are the eigenstates of 3^ m of the nonzero eigenvalues, 

m—1 

= ^(-in^M 172 n + ( m - r ) 2 ) $ t' y^% m = o . (4.24) 

The appearance of imaginary eigenvalues in the spectrum of the integral 3^.m will be discussed 
later. The action of y^ m on the states of the continuous spectrum (14.61) is given by 

m—1 

3W*t = T(-l) m kH(k 2 + (m-r) 2 )<S>± k , (4.25) 

r=0 
m—1 

y hm *t k = ±(-iy n kH(k 2 + (m-r) 2 )^t k , (4.26) 

r=0 

i.e. this integral, unlike the A^ m , see (14.171) . distinguishes the waves coming from the left 
and from the right, and as Xi <m , detects a difference between the states with distinct values 
of the upper (sign) index. 

As we have seen, behind the existence of the integrals of motion X^ m and 3^,m is the fact 
that there are two different Darboux-Crum transformations, which intertwine the Hamilto- 
nians H^ m and Hf m . In the case of the discrete operators 1Z[ and TZ m (one can also consider 
the intertwining operators V, T and T±, see the discussion in Section [2]), their composition 
transforms into a symmetry operation for the Hamiltonians H^ m , 

[1W m ,H±J = 0. (4.27) 

For extended system (14.11) . this composition corresponds to the integral 7Z = TZ{IZ m = 
diag(Tl{R m ,Tl{R m ), 

[n,ni, m ] = o. (4.28) 
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The composition of the intertwining relations generated by Xf m and Y^ m , (14.121) and (14.201) 
respectively, yields 

0^1%.-^- i,m) ^-1,171 = Hi,m(Yi*mX l m ), (X^ m Y lm )H \ m = H lm (X^ m Y lm ) . (4.29) 

The intertwining relations transform therefore into commutation relations, and correspond- 
ing integral of motion appears for each Hamiltonian, analogously to (14.27]) . The resulting 
integrals are the differential operators of the order In + 1 with n = m + I, and these are 
nothing else as the Lax integrals A^ n+1 in (I2.12p . which we rename here as 

= ^2n+l = Yl%iXi >m = Xj^Y^, \%%ni ^Lm] = • (4.30) 

For the extended system, these integrals of motion can be joined to form a diagonal operator, 
Z^ m , which is generated by the anticommutator of the previous conserved quantities, 

Zl,m = i (\ Z~ ) = \ {Xl > m ' ^ ■ (4 ' 31) 

The origin of the Lax integrals in the present extended Hamiltonian from the the inter- 
twining operators Xf m and Yf^ is illustrated on Fig. H] 

It is straightforward to check from the above relations that Z^ m commutes with the 
Hamiltonian (14. 2 p and with the integrals X^ m and 3^, m , 

[Z hm Xi, m ] = 0, [Z hm , y hm ] = . (4.32) 

Its square produces a polynomial in T-Li im , 

m-l [J- 1/2] 

V z (ni, m ) = Ui, m n n (Hi, m + (m - r) 2 ) 2 (^, m + (l-s- 1/2) 2 ) 2 . (4.33) 

r=0 s=0 

whose roots are all the discrete doubly degenerate energies of the extended system. Note 
that the roots of the bound states are of degree two, while the zero energy state root has 
degree one. All the corresponding energy eigenstates are the zero modes of the integral 2>i )fn , 

= 0, Z hm Q% m = 0, (4.34) 

which also detects the direction of propagation of the waves of the scattering sector, 

m-l [i— 1/2] 

Zi, m $i k = -(-lr^JI J] (k 2 + (m-r) 2 )(k 2 + (l-s-l/2) 2 )^ k , (4.35) 

r=0 s=0 
m-l [1-1/2] 

Zi,r**± k = {-l) m+l k\[ H ((k 2 + (m-r) 2 )k 2 + (l-s-l/2) 2 )<S>t k . (4.36) 

Because of the "PT-odd nature of the operators Af m and Bf m , from which Xi m , y^ m and 
Zi }Tn are composed, all this triplet of the integrals is PT-even. Thus, instead to be Hermitian 
operators, all the conserved quantities commute with the VT operator, 

[U ltm , VT] = 0, [X hm , VT] = 0, [y lim , VT] = 0, [Z hmt VT]=0. (4.37) 
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Fi gure 4: The "intertwining paths" , which correspond to the action of the operators X 2 i 
and Y 2 \ between the mutually shifted, self-isospectral Hamiltonians (presented by circles 
and triangles), are shown for the case of the extended PT-symmetric system ^2,1 = "Hs 3. 
The effect produced by these intertwining operators on the integer- valued lattice is the same 
as that of the discrete operators 7Zi and TZ m , respectively. On the half-integer- valued lattice, 
the horizontal and vertical distances between the corresponding systems are interchanged, 
and in this sense, the action of the operators X 2l and Y 2 \ on the half-integer-valued lattice 
is dual to that on the integer- valued lattice. Starting from any point, the composition of 
horizontal and vertical paths to the equivalent point on the same lattice produces the Lax 
integrals Z 2l and the symmetry operator 1Zi1Z m . 



It is this property of the PT-symmetry that requires the presence of the imaginary unit as a 
multiplicative factor in the definition of 3^, m in Eq. (14.211) . In turn, the factor i in Eq. (14.211) 
emphasizes then the existence of the splitting of the discrete eigenstates into two different 
families, and reveals an additional specific feature of the whole supersymmetric configuration 
we have here. The simultaneous requirement of a common basis of eigenstates for all the 
integrals of motion in addition to the relation (14.371) fixes that only the integral y^ m has 
imaginary eigenvalues for the set of doublet states $^ in Eq. (14.41) (see Eq. (I4.24p and also 
the examples in the next section). As all the integrals are mutually commuting operators, 
the picture here is different from that in a usual quantum mechanics where Hermitian (self- 
adjoint) mutually commuting operators possess a common basis of eigenstates with real 
eigenvalues. It is instructive to look in more detail what happens here. Explicit form of the 
states in (14. 4p shows that they are not the eigenstates of the VT operator. Indeed, they 
satisfies the relation 

pr*t = ■ ( 4 - 3g ) 

Remembering that X^ m and Z[ m annihilate the set of doublet states (I4.4p . while the extended 
Hamiltonian ( 14. II) has an entire real spectrum [and the states (14.41) are its eigenstates], one 
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concludes that the PT-symmetry has a broken nature just for the integral y^ m [we remind 
parenthetically here that the eigenstates from the continuous part of the spectrum have real 
eigenvalues for 3^,m, see Eqs. ( I4.25P and (I4.26p ]. Taking into account independently only 
the Hamiltonian operator, one can find another basis where these states are simultaneously 
the eigenstates of the Hamiltonian and the VT operator; therefore, for the TL^m the VT- 
symmetry is unbroken. 

We have identified the nontrivial integrals of the extended system and discussed their 
properties. Now we consider the related nonlinear supersymmetric structure. The diagonal 
matrix er 3 is a trivial integral of motion for Hi jm . Nevertheless it allows us to double the set 
of the nontrivial integrals of motion since the multiplication of any of them by 03 gives a 
new, linear independent nontrivial matrix integral of motion. So, in this way we obtain the 
set of six linearly independent nontrivial matrix integrals of motion 

G£ = *,, m) QfX = a z Q?X, (4.39) 
= yi,m , S\% = a 3 S^ , (4.40) 
^ = Z, m , CV> = a*dg n . (4.41) 

Notice the absence of the imaginary factor % in the definition of the second anti-diagonal su- 
percharges in comparison with the usual SUSYQM approach with a Hermitian Hamiltonian. 
This guarantees that all the three new integrals are also "PT-symmetric operators. 

The square of the matrix integral 03 equals 1, and it can be identified as the grading 
operator, T = <t 3 . This grading operator classifies then the Hamiltonian TLi )m and integrals 
^fmi a = 1,2, as bosonic operators, while the integrals of the antidiagonal matrix form, 

Qim and are classified as fermionic operators. In correspondence with this, we get a 
nonlinear superalgebra with the following set of nontrivial (anti)-commutation relations: 

{Qit egj = (-l) a+1 2S ab P x , = (-iy +1 25 ab P y , {Qi%,S}%} = -25 ab C^, 

[Qtl O = (-l) a 2e ab Si%P x , £gj = {-lY2e ab QflP y , (4.42) 

where P% = Px(Hi,m) an d Py = Py(Hi,m) are the polynomials defined in ( I4.15p . (I4.23P and 
(I4.33p . respectively. Note that the integral commutes with all the other integrals and, 
so, plays here the role of the bosonic central charge. 

The choice of 03 as the grading operator is, however, not unique. Another possibility 
corresponds, for instance, to the choice T = VT + (or, T = VT ). Indeed, this opera- 
tor is a (nonlocal) integral of motion, whose square is equal to 1. Such a grading oper- 
ator classifies the integrals as bosonic integrals, while and £^ are classified as 
fermionic integrals. The corresponding superalgebraic relations can be computed then by 
making use of the relations described above. In this case we have, particularly, a relation 
{^l^, Cfln] = ZSabPziH^m)- This corresponds to the fact that each of the unextended VT- 
symmetric systems Hf m and Hf m is characterized by the bosonized supersymmetry, in which 
the "PT-symmetric integrals iZ^ m and iZf m , respectively, are treated as the T = "PT + -odd 
supercharges. 

We summarize the whole picture on which the tri-supersymmetric structure is based on 
Fig. El 
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Figure 5: The basic blocks of the tri-supersymmetric structure of the extended VT- 
symmetric Hamiltonian Hi, m . Horizontal arrows inside a table correspond to the duality 
transformation induced by the symmetry (l,m) — > (m + |,Z — |). 

5 Examples 

To illustrate different properties of the systems H^ m and the tri-supersymmetric structure 
°f Hi,mi nere we present some examples for specific values of / and m. Before doing this, 
we first note that according to the relation (12.141) . the simplest nontrivial case of H^ (x) = 
Hi (x + iir) reduces, up to rescaling, just to the displaced reflect ionless Poschl- Teller system 
with one bound state in the spectrum. The unique bound state corresponds to a resonance 
with a complex energy value in the spectrum of the Hermitian Hamiltonian with real Scarf 
II potential V(x) = sech 2 (a;) — sinh x sech 2 x, which is depicted on Fig. [2j So we will consider 
more rich cases of reflectionless "PT-symmetric systems with two and three bound states. 

5.1 Systems with two bound states 

Without loss of generality, the family of reflectionless potentials V^ m with two bound states 
in the discrete part of the spectrum can be presented by non-negative integer values of 
the parameters I and m subjected to the condition I + m — 2. The case (I = 0, m = 2) 
corresponds here to the well known Hermitian reflectionless Poschl- Teller potential Vq 2 = 
Vq 2 — — 6/ cosh 2 x. By virtue of (12.14[) . this potential shares, up to rescaling, the same 
spectrum as complexified Scarf II potentials with (I — l,m = 1). Then the remaining case 
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H. 



sinhx 



2.0 



provides a first nontrivial example which is not related to a Hermitian (reflectionless Poschl- 
Teller) counterpart by means of shifting and rescaling of the coordinate. The potentials 
V = V 20 are solutions of the nonlinear s-KdV2 equation 



V {5) - 10V"'V - 20V'V" + 30V 2 V - 10(V" - 3VV) + W = 



(5.2) 



where = d 5 V/dx 5 . Notice that in contrast with the complexified case V^, the Hermitian 
version of the potential, V^o, which is plotted on Fig. [2J is not a solution of the s-KdV 2 
equation. The Hamiltonians H^ fall into the class of systems studied in in the context of 
the Vl~— symmetric nonlinear integrable systems; the corresponding extended Hamiltonian 
%2,o appears as a particular case of the diagonalized squared Dirac equation for a free spin- 
1/2 field in de Sitter space, see Ref. [37]. 

The degeneracy of the spectrum of 7^2,0 is twice that for each system in (15. 1 1) , and we 
have two doublets of bound states and two zero energy states at the very bottom of the 
four-fold degenerate continuous part of the spectrum. The bound states correspond here to 
the complex energy resonances in the Hermitian version with real Scarf II potential, see Fig. 
2. Their eigenf unctions, 
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(5.3) 



satisfy equations 



(5.4) 



and correspond to wavefunctions (14.41) . 

These solutions can be obtained from the non-physical states of the free particle, by 
applicacion of the operators (13.111) or (13.121) with I = 2 and m = in the case of the 
subsystem H^q- For generic values of I and m, the Darboux-Crum transformations that 
map the free particle eigenstates into those for the lower Hamiltonian H[~ m are realized in 
correspondence with intertwining relations ( 15.5ft by means of the operators 



vL = -v(prjv = -nvrjr . 



(5.5) 
(5.6) 



Note that in correspondence with relation H t , 
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in), operators (15.51) are ob- 



tained equivalently from the intertwining operators T>^ m and 2?j~ m by the half-period shift. 
The non-physical states which are transformed into bound states by means of the Darboux- 
Crum transformations ( 13. 1 1ft . ( 13.121) . ( 15. 5ft and ( 15.6ft are 
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and 

0+1/4 = cosh (|) » Ci /4 = sinh (~J , (5.8) 
which obey the same Schrodinger equations as in (15.4)) . 

9 1 

#0 4>% /4t = -4^-9/4' H 0_ 1/4 = -4^-1/4 • ( 5 ' 9 ) 

One can choose solutions of different parity with respect to P in (15. 1\ and (15. 8p to obtain 
the bound states (15. 3p . Choosing the functions with positive P-parity we have, 

< - K ±S»>v = f ( J| > v <*■«» 

where P 2 = ^ 2,0^1, 0; ^2 ,0 = ^2,0^-1,0^0 j ^2,0 = ^0,-1-^-1,1 > ^2,0 = ^0-^-2,0-^1,0 ■ Ex- 
pressions ( 15. 10p and (15. lip with the non-physical states of P-negative parity, i.e. <t>Zgu an d 

0~ 1 ^ 4 , remain almost identical up to multiplicative constant factors. One can choose these 
factors pure imaginary to produce, for each entry of ( 15 .101) and ( 15.111) . a state of definite PT- 
parity by starting from </>^ 9 / 4 or 4"-iu- This can be understood by taking into account that 
the intertwining operators have no definite P-parity, but they have a definite PT-parity. In 
general case, while the T>f m and T>\ m are PT-even, the operators Vf m and T>\ m are PT-odd. 

Using the same procedure as with bound states, we can construct the eigenstates in the 
scattering sector ( 14. 6 p by applying the intertwining Darboux-Crum operators to the plane 
waves (13. 15[) of H , 

W 2 ,„*J t = fc 2 4>^ , H 2 ,^t t = * 2 4>? t . (5. 13) 

The eigenfunctions ( 15. 12)) for H^q and if 2 ~o may also be obtained by applying, instead, the 
operators P^" and f>\ , respectively, to the same plane wave eigenstates. In the case of the 
zero energy eigenstates we have 

T")~ \ / T) - \ Q / n~ 2i arctan sinh x 

2,0 1 _ ^2,0 \ o 1 e 



$0 ~ { ±^|,o ) 1 ~ V ±^2,0 J ^ ~ ~4 ±e 2 * arctansinh * J ' (5 ' 14) 

i. e. the operators P^o and T>\ should act on the non-physical zero energy solutions of Hq 
which are proportional to x. Another solution of zero energy of H , which is a constant, is 
annihilated by the Darboux-Crum operators, 

^ol = ^2,ol = 0. (5.15) 

The extended Hamiltonian "H 2 ,o possesses three basic conserved quantities in the form of 
the matrix differential operators. One of these integrals, X 2t o, is given by 

-^2,0 \ _ I ^2,0^1,0^0,0^-1,0 



X 20= l " "2,0 J — j ^2,0^1,0^0,0^-1,0 (5 16) 

^2,0 / V ^-1,0^0,0^1,0^2,0 
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The explicit form of the higher order differential operators is 



*£> = & ± 7^t x & + ( 6 =F ^ ^hx - | cosh 2 x) £ (5.17) 

1 flO cink ^ -L „ - fl ^cX 2 ^ _ Ql "\ d _ 13 / 1 Q -I— Q/i' c^nk ^ /1 -L- n ci-^U ^ 3 



cosh^~x 



(12 sinh x ± i [| cosh 2 x — 3] ) ^ — ^4 - (13 =F 3i sinh x) (1 ± i sinh x) 



This integral acts on the physical states of the Hamiltonian "H 2 as follows, 

Af 2 ,o$J, = 0, ^,0^ = 0, A$y,^ = ±^$± (5.18) 

^i fc = ±(fc 2 +^)(^ 2 +^i fe , ^, = ±(^ + ^(^ + 0^. (5.19) 

The operator A^o does not distinguish the waves coming from the left or the right, but 
recognizes the states which correspond to resonances in the Hermitian Scarf II potential 
spectrum, by annihilating all of them. Another integral of motion is also an anti-diagonal 
matrix differential operator, 

where 



x 2,0 u / V -"2,0 



Y?o = B to = 3- ± • (5-21) 

2,u 2,u coshx v ' 

The integral of motion 3^2,0 commutes with A^o and encodes the information to be comple- 
mentary to that provided by the latter. This can be seen from its action on the physical 

states, 

3? 1 

y 2 ^% = ±-*± , y 2 , $± = ±-*± , 3^ $f = 0, (5.22) 

3> 2 ,o$t fe = Tk$%, y2,o$t k = ±k*± k • (5.23) 

The remaining doublet states are annihilated by 3^2,0 , which in this case correspond to the 
states of the zero energy at the bottom of the continuous spectrum. The waves coming from 
the left or the right are recognized by it, and like (15.161) . the 3^2,0 detects also the upper index 
of eigenstates. Note that the bound eigenstates here that correspond to resonances in the 
Hermitian Scarf II systems spectra have pure imaginary eigenvalues of 3^2,o- 

Coherently with the properties of the displayed antidiagonal integrals, the diagonal oper- 
ator i?2,o = ^2,o3^2,o annihilates all the doublet states, and separates scattering states coming 
from different directions, 

JZ 2 ,o$± = 0, Z 2 , $± = 0, Z 2)0 <S>± = , (5.24) 
3M>*J fc = -* (V + \) {k 2 + |) *± fc , Z 2 , *t k = k ^k 2 + 1) (V + ^ . (5.25) 
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5.1.1 Systems with three bound states 

Reflectionless systems with three bound states are constrained to fullfill the relation l+rn = 3. 
Potentials ([H]) with (I = 0,m = 3) and (I = 2,m = 1) are related by Eq. (pTTIj) . When 
(/ = 3, m = 0), the systems if^ have three bound states, all of which correspond to resonance 
states in the Hermitian version. So, all the mentioned three bound states systems have 
corresponding analogs in more simple cases we have discussed above. A more non-trivial 
example with three bound states is given by the Hamiltonians 

„. d 2 7 sinhx 

Hf 2 = -— ±5^_^. (5.26) 

ete cosh x cosh x 

The Hermitian counterpart potential for H^~ 2 is shown on Fig. 2. Here the potentials V = V^ 2 
satisfy the s-KdV 3 equation 

V {7) - UV(V {5) - 20V" V) - 42V (4 V - 70V'" V" + 70V"'V 2 + 70V' 3 - 140 V'V 3 + 
-21(V (5) - 10V'"V - 20V'V" + 30V V) + 84(V" - QVV') - 64V' = . (5.27) 

The extended Hamiltonian Hip has six bound states, where four correspond to the bound 
state from the set (13.11) . 



arctan sinh x 



0,m 



cosh x 

A arctan sinh x 



V 



J 



( 



\ 



3 =pi arctan sinhx 

cosh 2 X 

„=pi arctan sinh x 



-4$ 



± 



cosh 2 x 



(2 — 3i sinh x) 
(2 + 3i sinh x] 



\ 



) 



(5.28) 



(5.29) 



cosh 2 x 

These solutions are analogs of the bound states for the Hermitian counterpart system, which, 
in addition, admits resonances for a complex value of energy. Those resonances correspond 
to two extra bound states solutions in the spectrum of the Hamiltonian %\ 2, 



/ g—^- arctan sinh a; \ 



cosh 1/2 



X 



±i- 



g 4^ arctan sinh a; 



(5.30) 



/ 



\ cosh 1//2 x 

All the bound states described above, $o,m> $i,n 
physical states of the free particle, 0^ 4 , and 
solutions of i? are 



'PI 



and 



cosh 2x, 



coshx, 



and $o,/ ; can be derived from the non- 
y A , respectively. Here, the unphysical 

= sinh 2x , (5.31) 



sinh x : 



(5.32) 

they have the same eigenvalues (of H ) as the bound states <3>om an d m A5.29[) . The 

mapping between the states is given by 



0,m 



2i 
15 



tv{ 2 



-4 1 



± 



±^!,2 



01!, 



15 



^2 



-1/4 
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where we use the definition of the operators (13.111) and (15. 5p . Similar expressions can be 
found by making use of the operators (13.121) and (15.61) ; it is worth to note, however, that some 
non-physical states of the free particle cannot be mapped properly because are annihilated, 

V^Z 4 = £> 1)2 0I 4 = i>i,*4>U = ^, 2 0-4 = , (5.33) 
P+ 2 Ci = V{,2<P-i = A^Ci = i>l,2<P-i = . (5.34) 

The situation is quite similar to the previous case ( I5.15p . In fact, in this case the constant 
state of the free particle is also annihilated by the operators X^l = ^i ?\ = 0- 

The wave functions of the continuum are obtained by the same method from the free 
plane waves, 

;± ( ^1,2 \ Jkx ^± ( ^1,2 \ -ikx 



*«=U#,r • *-'=U#J e • (5 - 35) 

and have energies E = k 2 . 

The anti-diagonal, mutually commuting basic integrals of motion for this case, X\i and 
3^1,2, have differential orders |Af 1:2 | = 2 and 1 3^1,2! = 5, and read 

X Xfl\-( A 1,2 A 2, 



and 



Xh2 ~ [ x; a 0' )-{At t0 At >2 '0 ' ) ' (5 - 36) 



V -if Y W \ - if 5 l,2- B l J l- B l J 0- B l,-l- B l,-2 ^ /c o 7 X 

The explicit form of differential operators that compose (I5.36P and (I5.37P are 

Xt 2 = -7^7 ± -^7— -T- ^-5- (44 + cosh 2 x ± 20z sinh x) (5.38) 

dar cosh a: ax 8 cosh x 

and 

= -£ ± -^-A fV- (sinh 2 x ± 2z sinh 2) (V- (3 sinh x =f 2z ± i sinh 2 x) 

J-i^ aa: J cosh x ax 4 cosh x V / ax 13 cosh J x V 1 / ax* 



H (4 cosh 4 x + 10 cosh 2 x — 30 =F 65z sinh x ± lOz sinh x) 

cosh x \ / 

(2 sinh x =F 5i ± 4« cosh 2 x) . 



cosh 5 : 

The action of the integrals on the doublets of the Hamiltonian is given by 

Xifl*£ m = T™*t m , X ^tm = ±|*tn. ^2^ = 0, AT lia $± = ±i*J , (5.39) 

45? 

3>i l2 *J TO = °> 3> li2 $± m = 0, y lfl 9^ = ±32*£, 3^1,2^ = 0. (5.40) 

Note that, again, the states $ i) "which correspond to resonances in the Hermitian counter- 
part systems, are annihilated by the integral X\ 2 and are characterized by pure imaginary 
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eigenvalues of the second integral 3^i,2- The scattering states (I5.35P are eigenstates of the 
operators Af 12 and 3^1,2 , 

= =F (k 2 + $J fc , *i,2<^ fe = T (fc 2 + , (5.41) 

yi,2®% k = Tk(k 2 + 1)(A; 2 + 4)$J fc , 3^1,2$^ = ±*(* 2 + m 2 + 4)^ fc • (5.42) 

Finally, the diagonal integral Z\p. = ^1,2^1,2 = ^1,2^1,2 annihilates the whole set of doublet 
states, 

£i >2 ^ m = 0, Z 1>2 ^ m = 0, £ li2 $± = 0, 2 1>2 $± = , (5.43) 
and recognizes, as the integral 3^1,2, the waves coming from the left or the right, 

Zi,2$t k = + (k 2 + l)(k 2 + 4)<&± fc , = —k (k 2 + ^ (k 2 + l)(fc 2 + 4)$± fc . 

(5.44) 



6 Discussion and outlook 

In this paper, by analyzing a two-parametric family of reflectionless PT-symmetric Hamil- 
tonians, we have revealed a new supersymmetric structure. The class of potentials studied 
here provides an instructive example of quantum mechanical systems with non-Hermitian 
Hamiltonians. In comparison with the Hermitian version of the Scarf II potential, the spec- 
trum of its complexified counterpart contains two series of singlet states discovered earlier 
within a framework of the group theoretical approach. Surprisingly, this characteristic is 
imprinted in a tri-supersymmetric structure that is based here on the specific properties of 
the family of potentials: their pure imaginary period and discrete symmetries of a reflection 
type in the indexes. Usually, the imaginary period in both Hermitian and non-Hermitian 
Hamiltonians does not play explicitly an important role at the level of the spectrum, or in 
supersymmetric aspects. Following the original idea of Dunne and Feinberg for the case of 
a usual SUSYQM with a linear Lie super algebraic structure and mutually shifted (on a real 
line) Hermitian Hamiltonians [16], we construct an extended PT-symmetric system com- 
posed by two Hamiltonians with self-isospectral potentials, but now displaced mutually in 
the half of the imaginary period. The obtained composed system has three basic non-trivial 
integrals of motion, which in the generic case are the higher order differential operators. The 
importance of the splitting of the discrete states becomes clear by analyzing these integrals. 
Two of the anti-diagonal, super charge-type integrals, Xi >m and 3^,m, annihilate separately 
the two different sets of doublets of the extended system. These mutually commuting inte- 
grals generate a third, diagonal integral, Zi )fn , that implies the reflectionless property of the 
Hamiltonian: it appears as the Lax integral, which together with the Hamiltonian forms the 
Lax pair. The PT-operator emerges naturally as a valuable symmetry for the integrals of 
motion in view of the fact that all of them appear as PT-symmetric operators, in the same 
way as the Hamiltonian. Nevertheless, the odd order integral 3^,m reveals a finite number 
of pairs of complex conjugate eigenvalues when acts on the bound states which correspond 
to resonances with complex conjugate energy values in the Hermitian potential counterpart. 
We can say therefore that for this integral the PT-symmetry is spontaneously broken. 
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The "PT-symmetry is composed of the space inversion, V, and the time reversal, T, 
operators, which play a role of the intertwiners between the mutually displaced components 
ffj + m (i) and (x) = Hf m {x + in) of the extended Hamiltonian. In addition to them and 
the half-period displacement operators, there are other discrete intertwiners, the products 
of which produce discrete symmetries of the extended system. The commuting operators 
IZi : (I, m) — > (— l,m) and TZ m : (I, m) — > (/, —m — 1) produce, particularly, the same effect 
as the differential intertwiners X^ m and Y^, from which the supercharges X^ m and 3^ m are 
composed. In addition to a usual choice for the Z 2 -grading operator r = 03, other choices 
are also possible. The product of the V and of the operator of the displacement for the half 
of the imaginary period is one of them, which also happens to be the grading operator for 
the hidden, nonlinear bosonized supersymmetries of the subsystems H^ m and Hf m , where 
the Lax operators Zf m and Zf m are identified as the odd supercharges. 

On the other hand, the mentioned two sets of the discrete eigenstates can be related 
between themselves by means of another discrete symmetry of the Hamiltonian, which in- 
terchanges the integer-valued lattice of the parameters I and m with the half-integer-valued 
lattice. It is only for such, integer or half-integer, values of the parameters the complexi- 
fied Scarf II potentials are reflectionless. The indicated symmetry operation intertwines the 
operators Af m and Bf m , which are the building blocks for the intertwiners Xf m and Y^, 
respectively. As a consequence, the role of the integrals X^ m and y^ m is dually interchanged 
by those specific discrete symmetries. In contrast with the rest of the discrete symmetries, 
these duality generators have no analog in a form of differential operators. 

The supersymmetric structure presented here displays several similarities with the tri- 
supersymmetric structure in self-isospectral Hermitian finite-gap systems with elliptic po- 
tentials studied in [20]. In that class of Hermitian systems, two distinct finite-dimensional 
representations of si(2,R) are realized on periodic and antiperiodic band-edge states; like 
here three basic integrals of motion are present in the extended system, and different choices 
for the grading operators are also possible. The main difference with the present struc- 
ture is that the corresponding finite-gap elliptic systems are doubly periodic, in addition 
to the imaginary period the corresponding systems have also a real period, and there the 
self-isospectral systems are shifted for the half of their real period. The corresponding tri- 
supersymmetric systems studied in [20] are described by the associated Lame potentials, 
which constitute a subclass of the Darboux-Treibich-Verdier family f)2.2ip . It is interesting 
therefore to investigate the question of existence of tri-supersymmetric structure for such a 
class of doubly periodic 'PT-symmetric potentials, where the extended Hamiltonian would 
unify the self-isospectral partners with a mutual complex displacement. 

In the definition of a physically consistent, positively definite inner product for the sys- 
tems with non-Hermitian Hamiltonians, the existence of the operator C of the nature of a 
charge conjugation operator seems to be crucial [5]. An open question is the existence of 
such a kind of the operator for the complexified Scarf II potential. One can wonder then if 
the supersymmetric structure discussed here can be helpful in this sense, specifically, if the 
C can be expressed in terms of, or related to the non-trivial integrals of motion. 

Particular cases of the potential with I e Z and m = (and equivalent cases obtained 
by symmetry transformations of indexes) discussed here appear in quantum field theory in 
curved space-times [37J. A natural question is if a general case of the complex reflectionless 
potential plays any role in some related problems, and if the the revealed supersymmetric 
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structure could give some insight to these theories. 
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